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In this paper we investigate a simple parameterization scheme of the quintessence model 
given by Wetterich (2004). The crucial parameter of this model is the bending parameter 
b, which is related to the amount of dark energy in the early universe. Using the linear 
perturbation and the spherical infall approximations, we investigate the evolution of 
matter density perturbations in the variable dark energy model, and obtain an analytical 
expression for the growth index /. We show that increasing b leads to less growth of the 
density contrast S, and also decreases the growth index. Giving a fitting formula for the 
growth index at the present time we verify that the approximation relation / ~ also 
holds in this model. To compare predictions of the model with observations, we use the 
Supernovae type la (SNIa) Gold Sample and the parameters of the large scale structure 
determined by the 2-degree Field Galaxy Redshift Survey {2dFGRS). The best fit values 
for the model parameters by marginalizing on the remained ones, are — 0*21^q q^, 

wo = —2.05^2 Q5 ^^'^ ^ — 4.05^2 25 ^'^ confidence level. As a final test we calculate 
the age of universe for different choices of the free parameters in this model and compare 
it with the age of old stars and some high redshift objects. Then we show that the 
predictions of this variable dark energy model are consistent with the age observation of 
old star and can solve the "age crisis" problem. 

Keywords: theory dark energy; structure formation; observational constraints. 
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1. Introduction 

Recent observations of type la Supernovac (SNIa) at low and medium redshifts and 
the cosmic microwave background (CMB) anisotropics strongly indicate that the 
total matter-energy density of the universe is now dominated by dark energy i'2,3,4^ 
The origin and nature of this dark energy term remains unknown. The most obvious 
theoretical candidate of dark energy is the cosmological constant A which has the 
equation of state w = —1 Since the cosmological constant is a physical link 
between micro and macro scales, two main questions arise. First, the fine-tuning 
problem asks why the dark energy density today is so small compared to typical 
particle scales. The dark energy density is of order 10~^^ GeV^, which appears to 
require the introduction of new mass scale 14 or so orders of magnitude smaller 
than the electro-weak scale. The second difficulty, the cosmic coincidence problem, 
states, Since the energy densities of dark energy and dark matter scale so differ- 
ently during the expansion of universe why are they nearly equal today? To get this 
coincidence, it appears that their ratio must be set to a specific, infinitesimal value 
in the very early universe. As an attempt to solve these problems, one can consider 
a scalar field with a potential and kinetic term which has an equation of state to 
behave as vacuum energy. The energy density of this field, called quintessence in the 
cosmological attractor solution follows the energy dcmsity of matter and radiation, 
but remains negligible until recent epoch In these models the standard cosmolog- 
ical constant A-term is replaced by a dynamical, time-dependent component. Many 
different variable dark energy models have been proposed in the literature, e.g., Wet- 
terich (1988), Ratra & Peebles (1988) '^'^'^ For a more complete list of references and 
a review of this topic see This is still not a satisfactory physical explanation for 
the observed values of the cosmological constant term. The quintessence models do 
not offer a fundamental explanation, instead they are a phenomenological approach 
to express our inability of understanding the nature of the cosmological constant in 
terms of a variable scalar field. 

Here we examine a generic parameterization of quintessence models given by 
Wetterich (2004) in the recent paper some observational constraints have 
been investigated for this model and the value of present equation of state was 
assumed wq > — 1 but in this paper we let this parameter to hewo <0 and give the 
best fitting values for wa, present matter density and bending parameter. Also in 
the constraints related to the background evolution have been used. Our motivation 
for the present work is to see the possible observational effects I5'i6,i7 variable 
dark energy on the growth of the large scale structure of the universe. To put the 
rigorous constraints on the parameters of variable dark energy model we use the 
luminosity distance of Supernovae type la (SNIa) of the Gold Sample ^ and large 
scale structure parameters determined by the 2dFGRS team. The equation of state 
used for this quintessence model is: 



w{z;b,wo) 



[l + 61n(l + ^)] 



2 ' 



(1) 
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where wq is the state parameter at the present time and b is the bending parameter, 
which expresses the change in the equation of state of dark energy with redshift 
and is related to the amount of dark energy in the early universe. According to the 
theory of structure formation, the evolution of structure in the linear and non-linear 
regimes depends strongly on the background dynamics of universe. Dark energy as a 
crucial element of cosmic fluid aff'ects the dynamics of the universe and consequently 
changes the growing rate of structure. By increasing the parameter b in the model, 
the universe enters earlier in the phase of dark energy domination and the faster 
dilution of matter suppresses the formation of further structure. We assume that 
variable dark energy behaves as a smooth component, so that in our analysis the 
structure formation is only due to matter condensation, while the variable dark 
energy alters only the background cosmic dynamical evolution. We mostly work 
with 0(ot = 1 which is based on the results of CMB experiments 

This paper is organized as follows. In section 2 we discuss the linear perturbation 
theory of cosmological structure formation applied to the case of a variable dark 
energy. Considering the parameters of the model, we compute the evolution of the 
matter density contrast and the evolution of the growth index. We apply the spher- 
ical approximation with different initial conditions for overdense and underdense 
regions and compare these results with the linear approximation in Section 3. In 
section 4, we constrain the parameters of model by using high redshift Supernovae 
type la (SNIa) of the Gold sample and the parameters of Large Scale Structure 
(LSS) determined by 2-degree Field Galaxy Redshift Survey (2dFGRS). As a con- 
sistency test we look at the predicted age of universe in the model by comparing 
with the age of high redshift objects. The conclusions are given in Section 5. 



2. Linear Newtonian structure formation with variable dark 

energy 

The dynamics of the universe is driven by the Friedmann's equations as: 

-) +^ = — (p™ + Pa), (2) 

Pra + = 0, (3) 

a 

PA + 3^PA[l + w'(a;6,wo)] =0, (4) 

where the dot denotes the derivative with respect to time, p.^ and p\ stand for 
the pressureless matter and quintessence component, respectively. Also w{(i\ b, wq) 
is given by Eq. (1). As usual. A: = 0, 1,-1 indicates a flat, closed and open spacial 
section. Dynamic equations for each fluid (equations (3) and (4)) are: 

Pm = p'La-' , p. = p°a-^(^+'^(«^'"-°)), (5) 

where ■w{z; b, wo) = wo/[l + b ln(l + z)] . Fig. 1 shows the behavior of •w{z; b, wq) as a 
function of redshift. We notice that most variation in w takes place at low redshifts 
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and that larger values of b correspond to a more dust-like behavior of dark energy at 
high redshifts. The redshift of equality of matter and dark energy strongly depends 
on the choice of parameter 6, as shown in Fig. 2, a universe with a large b enters 
the dark energy dominated regime earlier. 

The Hubble parameter for a universe composed of dark matter and dark energy 
is given by: 

H^{z- b, wo) = H^prail + zf + n^iz; b, wo) - [ntot - 1)(1 + z)% (6) 
with the definitions: 

"~ 3 

nA(z;b,wo) = ^^(l + z)3[i+-(^;''-o)]^ (8) 
k 

^tot = 1 - — 2 = + J^A, (9) 

where fix and Q„i are the ratio of the dark energy density and dark matter density 
to the critical energy density at the present time, respectively. 



(7) 
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Fig. 1. tD(2; b, wq) in terms of redshift for various bending parameters. Here we choose wq = — 1. 
The slope of the graphs are more sensitive to the bending parameter at low redshifts rather than 
high redshifts. 
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Using the continuity and Poisson equations in the expanding FRW universe, the 
evolution of density contrast, 5 = Sp/p in the hnear approximation is given by i^-^O; 

S + 2-d~ (vs'^V^ + ^nGp,n)6 = 0. (10) 
a 

In Eq. (10) the dark energy enters through its influence on the expansion rate 
H{a; b, wq). The vahdity of this hnear Newtonian approach is restricted to pertur- 
bations on the subhorizon scales but large enough where structure formation is still 
in its linear phase i^-^o. "We also assume that the sound horizon of dark energy is 
much larger than the wavelength of the perturbations, so we do not need to consider 
the clustering of dark energy. If the perturbation is larger than the Jeans length, 
A,7 ~ n^/^Vs/ \/Gprm then Eq. (10) for cold dark matter (CDM) density contrast 
reduces to: 

J + 2-(^-47rGp„(5 = 0, (11) 
a 

The equation that describes the evolution of density contrast with respect to the 
scale factor is: 

^ d6 
da? da 

The numerical solution of Eq. (12) in a Friedmann universe with variable dark 
energy evolving as given in Eq. (6), is shown in Fig. 3. In the matter dominated 
era, the density contrast 5 grows linearly with the scale factor, while we have a 
deviation from a linear behavior, when dark energy begins to dominate. For larger 



a 



a 2H{a;b,'Wo) 



a 



3HI 
2d^a^ 



a,nS = 0. 



(12) 




Fig. 2. Redshift of equality of dark energy and dark matter density as a function of the bending 
parameter b, for a fiat universe (Solid line) with the Qrn = 0.3 and Qx = 0.7, an open universe 
(dashed line) with flm = 0.2 and Qx = 0.5 and a closed universe (dot-dashed line) with f!m = 1-0 
and Qx = 2.0. Increasing the bending parameter results in a higher redshift of equality. The present 
state parameter is set to wq = —1.0. 
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b (see Fig. 2) the domination of dark energy occurs earlier and the friction term in 
Eq. (11) increase so the growth of density contrast decreases. The relative difference 
between the present value of the ACDM and the dark energy models is about 2.5% 
for 6 = 0.1 and 100% for 6 = 1. 

In the linear perturbation theory the peculiar velocity field v is determined by 
the density contrast ^^'^^ : 



v(x) = Ho 



J_ 

in 



siy) 



X y 



(13) 



d\nS 



(14) 



where the growth index / is defined as: 

d In a 

We use the density contrast S to compute the growth index of structure /, which 
is an important quantity for the interpretation of peculiar velocities of galaxies, as 
discussed in for the Newtonian and ^■^ for the relativistic regime of structure 
formation. For understanding the physical meaning of the growth index it is helpful 
to divide the second term of Eq. (11) (friction) by the third (Poisson), which shows 
that: / cx 2HS/4:nGpmS. 

According to equations (12) and (14), the evolution of the growth index is given 

by: 



df 
din ( 



'Hi 



+ 2a3 " 



+ riA(a; b, Wq){1 + 3w{a; b, wq))] 



(15) 
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Fig. 3. Evolution of density contrast as a function of scale factor for different bending parameters 
in the flat universe with Cm = 0.3, Vl\ = 0.7 and uiq = —1.0. 
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In Fig. 4 we plot the numerical solution of Eq. (15) for the growth index as a function 
of redshift for 6 = 0, 6 = 1.0 and b = 4.0, for various amounts of dark matter and 
dark energy at the present epoch. The growth index / has been calculated for models 
with a cosmological constant earlier and the case of 6 = corresponds to the ACDM 
model with non- variable cosmological constant, where our results are identical with 
the earlier results given by There a fitting formula was given for the dependence 
of / on rim for the case of a flat universe. The dependence on the cosmological 
constant was found to be very weak. Here we extend this analysis to models with a 
variable dark energy, including the bending parameter b and Wq. In agreement with 
the tendency found earlier, we also find that in our case the dependence on b is not 
significant at the present time. This fitting formula is as follows: 

nx,b,wo,z = 0)!=^ f^^" + [-0.030 + 0.019 exp(-6) - O.OSSwq] 

Q[1.047-0.1456°'^^ + (-0.104-0.306^-^)«;g] _ 

^ 100 ■ ^ 

At the present time the growth index shows mostly a dependence on the matter 
density of the universe flm. But the time evolution of the growth index depends 
strongly on the choice of the parameters and b. Fig. 5 shows the time evolution 
of the growth index / for different pairs of (Clrn,^\), where the sensitivity to the 
bending parameter is examined. The main effect is that increasing the bending 
parameter decreases the growth index for all cases of open, closed and flat universe 
and results in a lower abundance of structure at the present time. In Fig. 6 the 
result from numerical solution of Eq. (15) and from fitting formula, Eq. (16), are 
illustrated. This shows that the growth index is approximated by / ^ 

In the case of a closed universe, the growth index rises to a maximum and 
descends afterward. This temporal increasing of the growth index is due to the 
dynamical effect of dark energy that changes the sign of the acceleration from 
negative to positive and the slope of deceleration parameter q = —d/aH^. During 
the transition phase of the acceleration, the universe reaches its lowest expansion 
rate and the structures have an opportunity to grow almost exponentially. Increasing 
the bending parameter suppresses the bump of the growth index (sec Fig. 5), because 
the dominance of the dark energy occurs at earlier times and increases the friction 
term proportional to the expansion rate. 



3. Variable Dark energy and nonlinear structure formation in the 
spherical approximation 

A simple approximation for calculating the formation of cosmological objects is the 
spherical infall model 23,24,25 ^j^jg j^odel we consider an overdense (or an under- 
dense) spherical region, i.e. a positive density perturbation 5 > (or 5 < 0) with 
spherical symmetry. For large enough positive perturbations, this region becomes 
gravitationally unstable and grows to a bound structure seen in the universe today. 
In such overdense regions the gravity of the perturbations is able to stop the expan- 
sion, turn it around and finally make the particles collapse. This highly nonlinear 
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process is usually modeled with a number of simplifying assumptions Here we 
also apply this approximation to the case of a negative density perturbation 5 < Q 
and take into account a non-constant dark energy term. For 5 < Q only the ini- 
tial conditions are different and the dynamical growth of the spherical underdense 
region is governed by the same equation. Since in the underdense regions the prob- 
ability of forming galaxies, in particular luminous galaxies is lower, these regions 
correspond to voids in the observed galaxy distribution. 

We begin with the Friedmann equation with a quintessence term as given in Eq. 
(6). Until an initial time ti we assume the spherical region Ri to grow with the same 
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Fig. 4. Growth index versus redshift for pairs of cosmological density parameters Cm and Vl\ 
The panels refer to different values of the bending parameter b and uiq = —1.0. Left panel: b = 0.0 
corresponds to the case of ACDM model. Middle: b = 1.0. Right: b = 4.0 as an extreme case of 
the variable dark energy model. 
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Fig. 5. The effect of bending parameter b on the evolution of the growth index as a function of 
redshift in case of different pairs of and Left panel: A fiat universe with f2m = 0.3 and 
= 0-7. Middle: An open universe with = 0.1 and Q,x = 0.5. Right panel: A closed universe 
with Qjn = 1-0 and Ux = 2.0. 
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rate as the background universe: 

2 




-3[l+w{ai\h,wo)\ 



(17) 



where is the initial scale factor of the background. At the time U we introduce 
an initial density perturbation Si in the matter density = (1 + di)QrnO'j^^ ■ As 
long as there is no shell crossing, the matter content inside a shell is constant. As in 
Section 2, we assume the effect of dark energy to be on much larger scales than the 
forming structure. This means that the dynamics of the universe is changed by dark 
energy but at the scale of structures, such as galaxies and clusters of galaxies the 
pressure of dark energy is negligible and the energy is conserved. Using the energy 
conservation for the spherical region with a radius R and a mass M yields: 

GM' 



= 2 

the potential energy at every time is: 
GM 



E 



R 



R 



using M — Mi, we get: 



(18) 



(19) 



(20) 



therefore by using equations (17), (19) and (20) we obtain for the size evolution of 
the spherical region R in the Eq. (18): 



0.55 1 I I I 1 1 I I I 1 1 I I I I 1 1 I I I 1 1 I I 1 1 I I I 1 1 I I I 1 1 I I I I. 
0.54 - 

0.53 - - 



^0.52 - 
O 




Fig. 6. The growth index / at the present time as a function of the bending parameter b for a 
fiat universe with = 0.3 and Q,\ = 0.7. The solid hne is given by the numerical solution of Eq. 
(15) and dashed-line is given by fitting formula (Eq. (16)). 
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_ ^m(l + Sj) Rl n2--i[l+w{a;b,wo)] 

-^2„,<S,4-("*ot"l)( — (21) 

After the time ti the scale factor a evolves according to the Friedmann equation 
and the spherical region changes its size with Eq. (21). For the numerical solution of 
this equation we replace R by dR/ da. Fig. 7 shows the evolution of i? as a function 
of scale factor a. In agreement with the results obtained by the linear perturbation 
theory (Fig. 3), Section 2, we also see a suppression of the growth of overdense and 
underdense regions if we increase h. Although the effect on the dynamical growth of 
spherical underdense regions seems to be very small, its tendency is discovered by 
observational data. It has been shown for example in ^6:27^ thaX large voids have a 
lower abundance in the cosmological iV-body simulations than in the observations 
of the large scale in redshift surveys. The effect of a variable dark energy would 
slightly increase this discrepancy. 



4. Observational constraints 

In this section we examine the constraints on the free parameters of the model from 
the observational data of the large scale structure (LSS) and SNIa. Observational 
constraints from the CMB and the amplitude of baryonic oscillations peak from the 
luminous red galaxies in Sloan Digital Survey Sky (SDSS) can be found in 




Fig. 7. Evolution of radius of a spherical overdense or underdense perturbation as a function 
of scale factor for different bending parameters. The initial conditions are Zi = 1100 and <5i as 
denoted in the diagram. 
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4.1. Constraints from Supernovae Type la 

The SNIa experiments provide the most important indication for the existence of 
dark energy in the standard model of cosmology. Since 1995 two teams, the High-Z 
Supernovae Search, and the Supernovae Cosmofogy Project have discovered several 
SNIa candidates at high redshifts 28,29 4 ^j^g discovery of 16 SNIa the Hubble 
Space Telescope were announced containing some of the most distant {z > 1.25) 
SNIa known to date. Based on this data a uniform Gold sample of high and low red- 
shift SNIa was constructed, 4,30,3i j-^ ^j-^jg subsection we compare the predictions of 
the dark energy model with the SNIa Gold sample. The observations of Supernovae 
measure essentially the apparent magnitude m including reddening, K correction 
etc, which is related to the (dimensionless) luminosity distance, of a an object 
at redshift 2;, for a spatially flat universe by: 

m = M + 5logDLiz;b,wo), (22) 

where for arbitrary spacial geometry 

Dl{z; b, wo) = ^"^f^ sinn -— 1 -dc) , (23) 

./iShTl V Jo H{C;b,Wo) J 

where ftK = 1 — ~ ^^a, and "sinn" is "sinh" for fc > (closed Universe) and 
"sin" for A; < (open Universe) For k = 0, Eq. (23) reduces to iJo(l + z) times 
the integral. Also 

^ = M + 51og,o(3^^)+25, (24) 

where M is the absolute magnitude. The distance modulus, /u, is defined as: 

/X = m - M = 5 logio Dl{z; b, wo) + 5 logio (p^^) + 25, (25) 

In order to compare the theoretical results with the observational data, we must 
compute the distance modulus, as given by Eq. (25). The first step in this sense is 
to compute the quality of the fitting through the least squared fitting quantity 
defined by: 

2 _ [Mobs (zj)- fJ-th {Zj ■,Clm,wo,b,h)]'^ 

i ' 

where CTj is the observational uncertainty in the distance modulus. To constrain the 
parameters of model, we use the Likelihood statistical analysis. The method and 
its motivation are described in detail in ref. . In the absence of prior constraint, 
the probability of the set of distance modulus jj, conditional on the values of a set 
of model parameters is given by a product of Gaussian functions: 

[liobs{zi) - ^ith{zi; 



pinthizz] {I})) = TT —7=^ exp 
i V2^ 



2a? 



(27) 



February 5, 2008 7:15 WSPC/INSTRUCTION FILE 

structure ' formation' 23july 



12 S Arbabi Bidgoli, M S Movahed and S Rahvar 



where {/} = {^rmWo,h,h}. This probabiUty distribution must be normalized. Evi- 
dently, when, for a set of values of the parameters, the is minimum the proba- 
bility is maximum. We find the minimum value of y^in/^d.o.f = 1-131 correspond- 
ing to the best fit values for the model parameters h = 0.66, wq = — l-SOtg 29, 
Q,m = O.0llo;oi and b = 6.00lg;Qo ^'^ level of confidence. As a special case we 
fix the value of the state parameter to wo = —1.0 and obtain the best fit values 
for the model parameters as: h = 0.65, Clm = O.OIIq q^ and b = 2.08lQ'gg at la 
confidence level with XmirJ-^d.o.f = 1.133. The values of cosmological parameters 
from fitting the dark energy model are different from those in the ACDM model, i.e. 
the value of h = 0.7ll^;^3 and f2„ = 0.27^^-^4 are sUghtly smaller in the WMAP 
results for ACDM model ^4,35 j-f ^^^^ Hubble parameter Hq = 71.0±7.0 from 

the HST-Key project as a prior parameter, the best fit values are: flm = O.OOIq'qo 
and b = 0.79l^:J^ at IfT confidence level with Xmin/^d.o.f = 1.552. The values of 
model parameters with Icr and 2<t confidence levels are summarized in Table 1 and 
2. Figures 8, 9 and 10 show marginalized relative probability density functions for 
cosmological parameters. 



4.2. Combined Constraints from Structure formation and SNIa 

The above analysis shows that the SNIa data alone does not sufficiently constrain 

the variable dark energy model. Furthermore the fit of SNIa data is sensitive to the 
parameter Hq. Hence, it is very important to find other observational quantities 
independent of Hq as a complement to the SNIa data. In the pervious sections 
we studied the effect of varying dark energy model on the evolution of large scale 
structures. Here we use the results from the 2dFGRS data and combine this with 
the SNIa data to put more rigorous constraint on the parameters of variable dark 
energy model. 

The 2dFGRS contains the observation of the position and the redshift of about 
220, 000 galaxies. On this basis the 2dFGRS team has determined a growth index /, 
which is the relevant parameter for the comparison with the predictions of the dark 
energy model Prom measuring the two-point correlation function, they report 

the redshift distortion parameter /3 = f /h = 0.49 ± 0.09 at z = 0.15, where h is the 
bias parameter describing the difference in the distribution of galaxies and mass. 
Verde et al. (2003) used the bispectrum of 2dFGRS galaxies ^^'^^ and obtained 
by^rde = 1.04 ± 0.11, from which, we obtain f{z = 0.15) = 0.51 ± 0.10. 

We perform a combined analysis of the SNIa and Large Scale Structure (LSS) 
to constrain of variable dark energy model through fitting: 

=XSNIa + XLSS, (28) 

where xiNia given by Eq. (26) for SNIa data and Xlss the contribution of LSS 
data using Eq. (15). The best fit values for the model parameters are: h = 0.66, 
Qrn = 0.2ltg;^^, b = 4.051^:^^ and wo = -2.05^^;^^ at la confidence level with the 
corresponding Xmin/^d.o.f = 1.131. The values of model parameters for la and 2a 
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confidence levels using SNIa+LSS observations arc given in Table 1. Fig. 8 shows 
the marginalized relative likelihood functions (upper panel) and joint confidence 
contours (lower panel) of the model parameters. 

Further, we restrict the analysis to a subset of the parameter space and fix the 
present state parameter at wq = —1.0. Again using the combined LSS and SNIa 
data, the best fit values for the model parameters are h = 0.64, fi^ = 0.20lo!o5 
and b = 0.80^°;go with X^„,„/A^d.o./ = 1-142 at la confidence level. The values of 
the model parameters for la and 2a confidence levels are reported in Table 2. Fig. 
9 shows the marginalized relative likelihood functions of matter density and the 
bending parameter. If we use Ho = 71.0 ± 7.0 from the HST-Key project as a prior 
with 1(7 measurement, we obtain the best fit values of the model parameters at 
flm = 0.101^:^^3, b = 0.20i2;2o at la confidence level with Xmin/^d.o.f = 1-618. The 
marginalized relative likelihood functions of the matter density and the bending 
parameter are shown in Fig. 10. 

Finally we determine the predicted age of universe considering the parameters 
of the model according to: 



and compare it with the age of old objects in the universe as a consistency test. In 
Table 1 we show that the age of universe from SNIa and the combined SNIa+LSS 
analysis are 13.45+^'^^Gyr and 13. 54^3:75 Gyr, respectively, which is in the range 
of the age of old stars 13^2 Gyr 39. 

Furthermore we use the observed age of three Old High Redshift Galaxies 
(OHRG) for comparison with the dark energy model, namely the LBDS 53W091, 
a 3.5-Gyr-old radio galaxy at ^ = 1.55 the LBDS 53W069, a 4.0-Gyr-old radio 
galaxy at z = 1.43 and a quasar, APM 08279 + 5255 at z = 3.91 with an age 
of i = 2.lto:?Gyr The latter one has once led to the "age crisis" in the ACDM 
model To quantify the age consistency test we introduce the expression r as: 



where t{z; b, wq) is the predicted age of universe from Eq. (29) and tots is an esti- 
mation for the age of the old object. In order to have a compatible predicted age 
of the universe, it is necessary to have r > 1. Table 3 shows the values of r for 
three mentioned OHRG. If we fix the parameters within the 1-a confidence level of 
the above combined observational constraints of SNIa and LSS, the predicted age 
of the universe is larger than the age of LBDS 53W069 and LBDS 53W091, while 
APM 08279 + 5255 at z = 3.91 is older. Only in the case that we fix = -1 and 
Ho = 71.0 ± 7.0 from HST-Key project, we obtain r > 1. 

5. Conclusions 




(29) 



T = 



t{z; b, Wo) 



(30) 



In this paper we have considered a simple model for variable dark energy and 
examined its influence on the cosmological structure formation. First we have used 
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the Newtonian linear perturbation theory to see the possible effects of a quintessence 
model, which is parameterized according to a prescription given by Wetterich (2004). 
The crucial parameter of this scheme is the bending parameter 5, which changes 
the behavior of dark energy toward a more dust-like behavior. Considering a flat, 
an open and a closed universe and an almost arbitrary initial value of the density 
contrast S, we find that by increasing the bending parameter, the growth of the S is 
decreased in comparison to the conventional ACDM model. In agreement with the 
results of our nonlinear spherical calculations, we find that increasing the bending 
parameter reduces the rate of growth of structure. This suppression in the growth 
of structure occurs in all open, flat and closed models for sufficiently large values of 
b. 

This effect can be explained by the fact that the dominance of the dark energy 
term occurs earlier with respect to ACDM, the larger the parameter b is. Because 





Fig. 8. The upper panel shows marginaUzed likelihood functions of cosmological parameters. The 
solid line was obtained by marginalizing over all nuisance parameters, using SNIa and dashed- 
line corresponds to SNIa+LSS. The intersections with the horizontal solid and dashed lines give 
the bounds for 68.3% and 95.4% confidence respectively. The lower panel shows joint confidence 
contours with Icr confidence level. 
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of the earlier domination of dark energy, the rate of the structure formation process 
is decreased. This finding is also supported by the second part of our calculations, 
where we numerically solve the equation of the collapse of a spherically symmetric 
density perturbation. We also give the fitting formula for the growth index, /, at 
the present time as a function of model parameters, the leading term of this formula 




Fig. 9. The marginalized likelihood functions of cosmological parameters. The solid line was 
obtained by marginalizing over all nuisance parameters, using SNIa and dashed-line corresponds 
to SNIa+LSS. The intersections with the horizontal solid and dashed lines give the bounds for 
68.3% and 95.4% confidence respectively. Here we assume wq = —1.0. 




Fig. 10. Marginalized likelihood functions of cosmological parameters. The solid line was obtained 
by marginalizing over all nuisance parameters, using SNIa+LSS+HST. HST prior is Hq = 71.0 ± 
7.0. The intersections with the horizontal solid and dashed lines give the bounds for 68.3% and 
95.4% confidence respectively. Here we assume uiq = —1.0. 
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is il,!^ which is in agreement with ACDM models 44,45,46 

In order to constrain the model parameters, we used the Gold sample SNIa data. 
The SNIa analysis resulted in a large degeneracy between the model parameters. 

For improving this, we perform a joint analysis of the SNIa and LSS data. This joint 
analysis of SNIa+LSS produces more reasonable results: h = 0.66, ilm = 0.2lto.o6i 
Wo = -2.05t^;[;^ and b = 4.05^1^1, at lo" confidence level with xLin/^d.o.f = 1-131. 
The age of universe from the best fit values of model parameters by using SNIa+LSS 
is 13.54t3"74Gyr, which is compatible with the age of old stars. We also chose three 
high redshift radio galaxies namely the LBDS 53W091, a 3.5-Gyr-old radio galaxy 
at z = 1.55, the LBDS 53W069, a 4.0-Gyr-old radio galaxy at 2; = 1.43 and a 
quasar, APM 08279 + 5255 at z = 3.91 with an age of t = 2.lt°;iGyr. The age of 
the two first objects is consistent with the age of universe using the best fit values 
of the model parameters, while the later one is older than the universe. The age of 
quasar APM 08279 + 5255, was also the reason of the "age crisis" for the ACDM 
model and it is compatible with the dark energy model only if we fix itjq = —1.0 and 
assume the Hubble parameter at the value given by HST-Key project as a prior in 
our calculations. 
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Table 1. Estimated variable dark energy model parameters (mean, 
68.3%, 95.4% C.L.), including the present matter density, the bend- 
ing parameter, the state parameter and the age of Universe. 
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Table 2. Estimated variable dark energy model pa- 
rameters (mean, 68.3%, 95.4% C.L.) for the fixed 
value of the state parameter uiq = — 1 , including the 
present matter density, the bending parameter and 
the age of Universe assuming wo = —1.0. The HST 
prior is set at 0.71 ± 0.07. 
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Table 3. Estimated relative age t of three high redshift objects using the best fit values 
of model parameters with Icr confidence level. 

Observation LBDS 53W069 LBDS 53W091 APM 08279 + 5255 
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